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SUMMARY 

The optimal control model (OCM) of the human pilot is applied to the study of air- 
craft handling qualities. Attention is focused primarily on longitudinal tasks. The 
modeling technique differs from previous applications of the OCM in that considerable 
effort is expended in simplifying the pilot/vehicle analysis. After briefly reviewing 
the OCM, a technique for modeling the pilot controlling higher order systems is intro- 
duced. Following this, a simple criterion for determining the susceptibility of an 
aircraft to pilot-induced oscillations (PIO) is formulated. Finally, a model-based metric 
for pilot rating prediction is discussed. The resulting modeling procedure provides a 
relatively simple, yet unified approach to the study of a variety of handling qualities 
problems . 


1. INTRODUCTION 

The advent of modern digital stability and control augmentation systems has created 
a renewed interest in the study of aircraft longitudinal handling qualities. This 
renewed interest is attributable to two factors: First, the higher order nature of the 

dynamics typically associated with digital control systems makes analytical prediction 
of handling qualities difficult. Contemporary handling qualities specifications (Ref. 1) 
are written assuming "classical” aircraft characteristics, e.g., in the longitudinal 
mode, the existence of distinct and dominant short-peTlod dynamics is assumed. With 
modern systems, the short-period characteristics may be dramatically altered by feedback 
and the higher order control system dynamics may dominate the vehicle handling qualities. 
Second, shortcomings in predictive techniques are made even more critical by the fact 
that severe handling qualities deficiencies often arise in practice which are directly 
attributable to the higher order nature of the digital control law implementation. An 
example of this is the ability of high frequency phase lags or time delays in the control 
system to sharply degrade aircraft handling qualities and to be a contributing factor to 
pilot-induced oscillations (RefT 2). 

In the research to be described, a pilot-modeling technique for handling qualities 
research, discussed in Ref. 3, is utilized and extended to cover higher order systems. 

The characteristics of over thirty aircraft configurations are analyzed, primarily in 
the longitudinal mode. Particular emphasis is placed upon those configurations where 
control system dynamics and time delays have been recognized as contributing factors to 
handling qualities deficiencies. The contribution of vehicle/control system dynamics 
to PIO tendencies is outlined, and a metric for pilot rating prediction is discussed. 


2. BACKGROUND 

The pilot-modeling technique as discussed in Ref. 3 forms the framework for the 
research described here. This technique utilizes the optimal control model of the human 
pilot and a novel method for the a priori selection of dominant OCM parameters (index 
of performance weighting coefficients and observation noise/signal ratios). A brief 
tutorial review of the procedure for selecting index of performance weighting coefficients 
is now presented. Consider the longitudinal tracking task of Fig. 1 in which the pilot 
is attempting to minimize pitch attitude deviations 6(t) in the presence of atmospheric 
disturbances. Ignore the dashed "internal attitude command" for the present. An accept- 
able index of performance for this task would be (Ref. 3) 


J ■ p(t)/e2 + 


where 6(t) is control rate. 

As discussed in Ref. 3, we assign an arbitrary maximum allowable deviation to the 
time rate of change of the error, $(t), and denote it 0^. Now an effective time constant 
T can be introduced to define maximum allowable deviations of the integral and derivatives 
of 0(t) as: 
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and 


0JJ ■ specified but arbitrary 


( 1 ) 


0m/t - e„/T" 


The Justification for using a single time constant to represent the ratio of the maximum 
value of a variable to that of its next highest derivative rests upon the system bandwidth 
implications which follow when Eq. (1) is used in implementing the OCM. We will also 
assign a maximum, allowable deviation to the time rate of change of the pilot’s control, 
i(t), and denote it Similar to Eq. (1) we write 



** to be selected ; 

: ( 2 ) 


The value of 6^ is not arbitrary, however, but is found using Eqs. (1) and (2) and the 
vehicle dynamics as follows: Let the pitch attitude dynamics of the aircraft be given by 


,n-l 


-n-2 


I (s) 


*n-2“ 


* * 




. + b-s + b^ 
1 o 


Then, as explained in Ref. 3, we write 


l/T°-^ ■. ^ |b^| ^ |b^|T . 

. |a^_2|/T“-3 ^ + | | + |a^|T) “ 


(3) 


(4) 


Thus, once T is known, and 0 m (and, if needed, 0y, etc.) can be determined immedi- 

ately. Choosing T involves selecting a domain of 1/T: 1/4t < 1/T < 4/x and then 
plotting J, the value of the OCM index of performance, vs 1/T. The operating point or 
"knee” of this curve determines T. The knee is defined as the point where 


?J „ *^It-t/4 ~ '^It"4t 

31og(l/T) “ ^^6 1^ f 4/ T ) log ( 1 / 4 T ) 


Here is a constant, nominally unity, which can be used to reflect manipulator 

characteristics, much like an efficiency factor; x is the pilot's time delay (nominally 
0.2 sec). J|t»x/ 4 value of the index of performance which results when T ** x/4. 

The ability of the OCM parameter selection technique to provide a pilot model which 
matches measured pilot describing functions, remnant power spectral densities and root 
mean square (RMS) performance measures was demonstrated in Ref. 3. In addition, the 
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modeling technique was shown capable of providing qualitative and quantitative handling 
qualities assessments. The method for selecting observation noise/signal ratios for 
the OCM is discussed in Ref. 3 and will not be dealt with here. 

Although Eq. (3) shows dynamics of arbitrary order, all the pitch attitude dynamics 
of Ref. 3 were of the form: 


e 

6 


Ke<s + l/T^) 

5 5“ 

sCs*^ + ^ 


( 6 ) 


When higher order dynamics are encountered, the method for selecting the operating point 
needs to be modified slightly. The large phase lags typically associated with the 
dynamics of vehicles with higher order dynamics need to be reflected in choosing the 
domain of 1/T to be used in Eq. (5). To accomplish this, a delay tq is defined as 
the delay which accrues when the vehicle dynamics of Eq. (3) are represented as 


0 

6 


Kq(s + l/\)e ^ 

' 2 
s(s + 2c„w^s + ) 

n n n 


(t) 


The parameters on the right hand side of Eq. (7) are found using a program to fit a 
linear transfer function model to the actual vehicle dynamics (Ref. 4). Equation (5) is 
modified by simply replacing t with t + tq. The resulting equation is interpreted 
graphically in Fig. 2. Calculating td and including it in Eq. (5) constitutes the 
extension of the methods of Ref. 3 to higher order systems. It is important to emphasize 
that the actual higher order vehicle dynamics are used in the modeling procedure; Eq. (7) 
is employed only to select td which, in turn, determines the domain of 1/T used in 
finding the index of performance weighting coefficients. 


3. APPLICATION TO AIRCRAFT HANDLING QUALITIES 

3.1 Pilot-Induced Oscillations 

Table I lists the aircraft configurations which have been analyzed in this study. 

The designations in the column labeled ’’Configuration" use notation found in the corre- 
sponding references. The first sixteen deal with high performance fighter-type aircraft 
in tracking or landing approach conditions and are taken from Refs. 2, 5, and 6. These 
configurations constitute the test cases for the majority of the assessments. The next 
four configurations are taken from Ref. 7 and represent pilot-in-the-loop simulations of 
a hovering helicopter. Configurations 21-25 are flight test results from Ref. 9 in 
which the Princeton University Variable Response Aircraft (VRA) was used to determine 
the effect of digital* sampling rates and time delays on longitudinal handling qualities. 
The vehicle dynamics appropriate for 105 kts airspeed were used in the modeling procedure. 
The pilot ratings used were average values obtained from altitude tracking and approach 
and landing tasks (Fig. 3 of Ref. 9). Finally, configurations 26-32 are taken from 
Ref. 10 where a moving-base simulator experiment on the NASA Ames Flight Simulator for 
Advanced Aircraft (FSAA) was described which investigated a wings level-turn control 
mode for air-to-ground weapons delivery. Note that unlike the previous twenty-five 
configurations, these involve lateral-directional aircraft handling qualities. The 
effective vehicle dynamics for the lateral gunsight aiming task were parameterized by 
a damping ratio Cn» undamped natural frequency Wq and a pure time delay tp 
(Ref. 10). The data for the so-called "fine" task were used. This task is explained in 
Ref. 10. 

As an example of the modeling results. Fig. 3 shows the longitudinal open-loop 
pilot/vehicle characteristics (YpYc) for three of the configurations used in Ref. 2. 

Here, the NASA Dryden F-8 digital fly-by-wire aircraft is considered with a rudimentary 
augmentation system ("Pitch Direct") and three transport time delays of 0.13 sec, 0.23 
sec, and 0.33 sec, respectively. The predicted effect of the time delays is apparent 
in the reduced open-loop crossover frequencies Wc* This open-loop characteristic 
obviously has a deleterious effect on the closed loop 0/9c transfer functions as shown 
in Fig. 4 [0/0c * YpYc/(l + YpYc)] . This transfer function is important in assessing 
PIO susceptibility. Although the task has been defined as pitch-attitude disturbance 
regulation, attitude commands 0^ Internally generated by the pilot would be employed 
in precise altitude regulation (dashed line in Fig. 1). Note in Fig. 4, that as td 
increases, |0/0cl and ^.0/0^ decrease at all frequencies. Perfect command following, 
of course, implies 0/0^ =• 1.0 at all frequencies. In Fig. 4, |0/0cl < 1*0 for all 
configurations when w < 3.0 rad/sec, and is particularly poor for the configuration 
with Td “ 0.33 sec. It can be readily shown that open-loop crossover frequencies less 
than 3-4 rad/sec will invariably result in poor closed-loop attitude command-following 
characteristics. The simplest and most direct way for the pilot to attempt to improve 
this closed loop command-following performance is to increase cjc by increasing his 
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static gain. If the pilot attempts this for the F-8 configuration with » 0.33 sec, 
a very lightly damped closed-loop oscillation occurs at w » 3.3 rad/sec (see Fig. 4). 

This is identical to the PIO frequency shown in Ref. 2 for this configuration. 

Similar results are also obtained for configurations from Ref. 5. Figure 5 compares 
a pair of open-loop transfer functions obtained using configurations "11" and "12" from 
Ref. 5 and applying the pilot-modeling technique discussed above. Once again, the 
dramatic difference in the crossover frequencies is apparent. The effects of the 

pilot's attempting to improve the performance of configuration "12" by increasing his 
static gain by 10 dB are shown in Fig. 6. Once again, a lightly damped oscillatory 
mode is seen to appear. The simulations of Ref. 5 were intended to provide performance 
comparisons for configurations which were flight tested and discussed in Ref. 8. The 
latter report included Pllot-Induced-Oscillatlon-Ratings (PIOR) obtained using the scale 
of Fig. 7. It is interesting to note that configuration "11" received an average PIOR 
of 1 indicating a very satisfactory vehicle whereas configuration "12" received a 
marginal average rating of 2.7 indicating a vehicle with definite PIO tendencies. These 
experimental results are seen to corroborate the analytical findings just discussed. 

Next, consider two configurations from Ref. 6 denoted as "4-1" and "6-1." Figure 8 
shows the YpYc plots for these configurations. Configuration "4-1" received a very 
satisfactory. PIOR of 1 whereas configuration "6-1" received a very poor PIOR of 4. 

Indeed, configuration "6-1" produced a PIO in flight test with a frequency of approxi- 
mately 3.75 rad/sec. Analytically increasing the pilot's static gain by 4.75 dB (the 
limit for closed-loop stability) in the modeling-results for this configuration produced 
a closed-loop oscillation at approximately 3.5 rad/sec. This 4.75 dB increase would 
increase from around 1.5 rad/sec to only around 2,5 rad/sec as compared to a value 

of 4.5 rad/sec for configuration "4-1." 

Figure 9 shows the predicted YpY^’s for a pair of configurations from Ref. 9. 

The task was longitudinal control in approach and landing using the Princeton VRA. 

The variable of interest here was the amount of effective delay in the control system. 

In the first, an effective delay of 0.055 sec was employed, while in the second, 0.355 sec 
was used. Again, note the striking difference in crossover frequencies in the predicted 
pilot/vehicle dynamics. In the first case, Wq *■ 3.4 rad/sec, while in the latter, 

» 0.55 rad/sec. Flight test of the first configuration showed no PIO tendencies, 
while those for the latter produced PIO's (Ref. 9). 

Finally, Fig. 10 shows the predicted YpY^'s for a pair of configurations from 
Ref. 10. In the first, the control system parameters were Cn ~ 1-4, Wn “ 2.0 rad/sec 
and Tp * 0.0 sec, while in the second, Cn “ 1-4, Ujj =* 15.0 rad/sec and » 0.49 sec. 

The u)c difference is again evident. Simulation results indicated that the configura- 
tion with delay was definitely PIO prone and the one without delay was not. It is 
interesting to point out that the configuration without delay still received an average 
Cooper-Harper pilot rating of 6.5, even though it was not PIO prone. Thus, poor pilot 
ratings, per se, are not a necessary condition for PIO susceptibility. 

In each of the cases above, we have made direct comparisons of vehicles which were 
found to be PIO prone with those which were not. This was done to emphasize the fact 
that the method proposed here is clearly discriminatory in predicting PIO susceptibility. 
The simple criterion for exonerating a vehicle from PIO tendencies requires that the 
predicted pilot/vehicle crossover frequencies associated with inner attitude-loops be 
greater than 3-4 rad/sec. 

3.2 Cooper-Harper Ratings 

Figure 11 is a plot of the Cooper-Harper ratings which the thirty-one configurations 
from Table I received in simulation or flight test vs the value of a proposed handling 
qualities metric defined as K4 C(t + td)/t1"*J. No ratings were reported in Ref. 10 for 
configuration 32 of Table I. Hence, only thirty-one data points are shown in Fig. 11. 

The Ki can be interpreted as a "calibration parameter" which, when multiplied by 
C(t + allows the reported pilot ratings from different tasks and data sources 

to coalesce as shown in Fig. 11. Note that we do not allow Ki to vary within the 
analysis of any particular task, regardless of configuration changes. Thus, the analysis 
of the six configurations from Ref. 2 used a single value of Ki (call it Ki). The 
analysis of the seven configurations from Refs. 5 and 8 used a single value of Ki (call 

it K2), etc. In all, six different Ki values (each one corresponding to the six 

different symbols in Fig. 11) were used to generate Fig. 11. With the exception of Ki , 
all the parameters of the metric are an intrinsic part of the modeling procedure, and, 
as such, involve no guesswork on the part of the analyst. In order to determine Ki, 
the analyst must have an actual pilot rating for one of the configurations tested for 
the task under study. If the analyst does not have such a rating available, Fig. 11 is 
still useful, since the curve is nearly linear from a pilot rating of about 2.0 to 10.0, 

a range which covers 80% of the Cooper-Harper scale. Thus, relative rating changes may 

be able to be predicted using the linear portion of the curve. Note that, with the 
exception of one data point (Config. 19 from Ref. 7), the scatter in the ratings in 
Fig. 11 is only about a pilot rating. 

The inclusion of the factor C(t + td)/!!** in the metric deserves a brief discussion. 
In previous research with the OCM, the value of J, alone, has been found to correlate 
well with pilot opinion rating (Ref. 11). In many of the configurations studied here, 
however (those with tj) > 0), the value of J was not acceptable as a metric. In 
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general, the ''predicted" opinion rating increments were smaller than those reported in 
experiment. There appears to be a reason for this based upon pilot tracking performance. 
Namely, when the task is disturbance regulation involving relatively low-bandwidth 
turbulence, large time delays are not necessarily a harbinger of dramatic deterioration 
in tracking performance. This is analytically verified by considering the RMS tracking 
scores for configurations 1 and 3 from Table I. Here, a 154% increase in time delay 
between configurations 1 and 3 involves a log uc regression of nearly a decade. 

However, the predicted RMS pitch attitude score Increases by only 36% and the predicted 
RMS control-rate score actually decreases. As we have attempted to point out here, 
however, the same cannot be said for discrete command following or abrupt maneuvers. 

In this case, wc regression can have a significant impact on the ability of the closed- 
loop pilot/vehicle system to follow abrupt, internally generated commands. It certainly 
is not unreasonable to postulate that such short-term response characteristics (in addi- 
tion to RMS characteristics) are reflected in pilot opinion rating. Indeed, recorded 
pilot comments support this notion (e.g.. Refs, 2 and 6). The inclusion of C(x + td)/!!** 
in the metric appears to account for the Influence of these delays on pilot opinion in 
a straightforward manner, employing an easily identifiable parameter (tjj). 


4. CONCLUSIONS 

The research summarized in this paper provides a unified approach to pilot/vehicle 
analysis, and in particular for: 

1) Modeling the pilot controlling higher order systems. 

2) Predicting the susceptibility of aircraft to longitudinal PIO's. 

3) Predicting pilot ratings for tasks when one configuration rating is known, or 
predicting relative rating changes between configurations. 

Although the majority of tasks studied dealt with longitudinal control, five lateral- 
directional configurations were successfully analyzed with no changes in the modeling 
technique. 
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Table I. Aircraft Configurations Analyzed 


No. 

Configuration 



Reference 

1 

F-8 "Pitch Direct" 

0.13 sec 

delay 


2 

2 


.23 



2 

3 


.33 



2 

4 

"ISAS" 

0.13 sec 

delay 


2 

5 


.23 



2 

6 


.33 



2 

7 

”2D" 




5,8 

8 

"5A" 




5,8 

9 

"8A" 




5,8 

10 

•’9” 




5,8 

11 

"10" 




5,8 

12 

"11" 




5,8 

13 

"12" 




5,8 

14 

"3-1" 




6 

15 

"4-1" 




6 

16 

"6-1" 




6 

17 

"PH-28" 




7 

18 

"PH-29" 




7 

19 

"PH-32" 




7 

20 

"PH-35" 




7 

21 

Princeton VRA 

0.055 sec delay 


9 

22 


.135 



9 

23 


.255 



9 

24 


.355 



9 

25 


,455 



9 


FSAA Wings-Level Turn 
(lateral-directional ) 


0} 

n 

( rad/sec ) 

(sec) 

10 

26 


1.4 

15.0 

0 

10 

27 


1.4 

2.0 

0 

10 

28 


2.0 

8.0 

0 

10 

29 


0.7 

6.0 

0 

10 

30 


0,5 

4.5 

0 

10 

31 


0.3 

4.5 

0 

10 

32 


1.4 

4.5 

0.49 

10 


INTERNAL ATTITUDE COMMAND 




Figure 1, A pitch attitude regulation 
task. 


Figure 2. Selecting an "effective time 
constant" T. 
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.1 1 10 


03 rad/ioc 

Figure 3. Pilot/vehicle dynamics for 
three configurations from 
Ref. 2. 


^ FOR To - 0.33 sac AND 5 dS INCREMENT IN Yp 



03 rad/sac 

Figure 4. Closed-loop characteristics 
for three configurations 
from Ref. 2. 



03 rad/tec 


Figure 5. Pilot/vehicle dynamics for 
two configurations from 
Ref. 5. 



03 rad/sac 


Figure 6. Closed-loop characteristics 
for a configuration from 
Ref. 5. 
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:OOPER-HARPER RATING 



Figure 9, Pilot/vehicle dynamics for 
two configurations from 
Ref. 9. 



Figure 10. Pilot/vehicle dynamics for 
two configurations from 
Ref. 10. . 


CONFIGURATION 





Figure 11. Cooper-Harper pilot ratings 
vs a proposed model-based 
metric. 
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